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Abstract 

In the framework of the /cy-factorization approach, we study the production of Higgs 
bosons associated with a heavy (beauty or top) quark pair at the CERN LHC colhder 
conditions. Our consideration is based mainly on the off-shell gluon-gluon fusion suprocess 
9*9* ~^ QQH . The corresponding matrix element squared have been calculated for the first 
time. We investigate the total and differential cross sections of hhH and ttH production 
taking into account also the non-negligible contribution from the qq QQH mechanism. 
In the numerical calculations we use the unintegrated gluon distributions obtained from the 
CCFM evolution equation. Our results are compared with the leading and next-to-leading 
order predictions of the coUinear factorization of QCD. 



PACS number(s): 12.38.-t, 12.38.Bx 



1 Introduction 

It is well known that the electroweak symmetry breaking in the Standard Model (SM) 
of elementary particle interactions is achieved via the Higgs mechanism. This mechanism 
is responsible for the generation of masses of the gauge {W^ and Z^) bosons as well as 
leptons and quarks via Yukawa couplings. In the minimal model there are a single complex 
Higgs doublet, where the Higgs boson H is the physical neutral Higgs scalar which is the only 
remaining part of this doublet after spontaneous symmetry breaking. In non-minimal models 
(such as Minimal Supersymmetric Standard Model, MSSM) there are additional charged and 
neutral scalar Higgs particles. At moment, the Higgs boson is the only missing, undiscovered 
component of modern particle physics, so that the search for the Higgs boson is of highest 
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priority for particle physics community. It takes important part at the Tevatron experiments 
and will be one of the main fields of study at the LHC collider [1]. The lower bound on 
the SM Higgs boson mass from direct searches at the LEP2 energy is mn > 114.4 GeV [2], 
while the recent global SM fits to electroweak precision data imply uih < 211 GeV [3]. The 
MSSM requires the existence of a scalar Higgs boson lighter than about 130 GeV, so that 
the possibility of Higgs discovery in the mass range near 115 - 130 GeV seems increasingly 
likely. 

The associated production of a Higgs boson with a heavy (beauty or top) quark pair can 
play a very important role at high energy hadron coUiders. At the LHC, the ttH production 
is an important search channel for Higgs masses below 130 GeV [4-6]. Although the expected 
cross section is rather small, the signature is quite distinctive. Moreover, analyzing the tiH 
production rate can provide information on the top-Higgs Yukawa coupling [6-9], assuming 
standard decay branching ratios [6] , before model independent precision measurements of this 
couphng are performed at e+e" colliders [10-12]. The Higgs boson production in association 
with two beauty quarks is the subject of intense theoretical investigations also [13-15]. In 
the SM, the coupling of the Higgs to a bb pair is suppressed by the small factor rrib/v, where 
V = (\/2Gf)~^^'^ = 246 GeV, implying that the bbH production rate is very small at both the 
Tevatron and the LHC energies. However, in the MSSM this coupling grows with the ratio 
of neutral Higgs boson vacuum expectation values, tan /3, and can be significantly enhanced 
over the SM coupling. Therefore it is one of the most important discovery channels for 
super symmetric Higgs particles at the LHC. 

From the theoretical point of view, the cross section of Higgs and associated heavy 
quark pair production at the Tevatron and the LHC is described by the gg — > QQH and 
qq — > QQH subprocesses (at the tree level). The leading-order (LO) QCD predictions [14, 
15] are plagued by considerable uncertainties due to the strong dependence on the renormal- 
ization and factorization scales, introduced by the QCD coupling and the parton (quark and 
gluon) densities. First estimates of radiative corrections were performed [16] in the so-called 
"effective Higgs approximation" (EHA). Recently the calculations of the 0{ag) inclusive 
cross section for the bbH and ttH production have been carried out at next-to-leading order 
(NLO) [17-20] of QCD. These calculations are based on the complete set of virtual and 
real 0{as) corrections to the parton level processes gg QQH and qq QQH, as well 
as the tree level process {q,q)g QQH + {q,q). The NLO cross sections are about 20% 
smaller and about 30% larger than the relevant LO cross sections at the Tevatron and LHC 
conditions, respectively. It was demonstrated [19, 20] that these high-order QCD corrections 
greatly reduce the renormalization and factorization scale dependence of LO results and thus 
stabilize the theoretical predictions. 

In the present paper we will study the Higgs and associated heavy (beauty and top) 
quark pair production using the so-called /cT-factorization QCD approach [21-24]. This 
approach is based on the famihar Balitsky-Fadin-Kuraev-Lipatov (BFKL) [25] or Ciafaloni- 
Catani-Fiorani-Marchesini (CCFM) [26] equations for the non-collinear gluon evolution in a 
proton. Detailed description of the /cT-factorization approach can be found, for example, in 
reviews [27-29]. Here we would like to only mention that the main part of high-order radiative 
QCD corrections is naturally included into the leading-order /c^^-factorization formalism. 

The /cT-factorization approach has been already applied [30-35] to study the inclusive 
Higgs production at the Tevatron and LHC conditions. First investigations [30-33] were 
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based on the amplitude for scalar Higgs boson production in the fusion of two off-shell 
gluons g*g* — H. The corresponding matrix elements have been derived first in [36] using 
the large rrit limit where the effective Lagrangian [37] for the Higgs boson coupling to gluons 
can be appliecj^]. The investigations [30-33] provocated further studies [34, 35] where the off- 
shell matrix elements of g*g* —>■ H subprocess have been calculated including finite masses 
of quarks in the triange loop. It was claimed [35] that the fc^-factorization approach give us 
the possibility to estimate the size of unknown coUinear high-order corrections. 

The starting point of present consideration is the off-shell amplitude of gluon-gluon fusion 
suprocess g*g* — » QQH. We evaluate the corresponding matrix elements squared for the 
first time and apply them for investigation of the bbH and ttH production rates at the LHC 
energy, a/s = 14 TeV. The quark- ant iquark annihilation mechanism, qq QQH, is expected 
to be significant only at relatively large x, and therefore we can safely take it into accout in 
the usual leading-order coUinear approximation of QCD. In the numerical calculations we 
will use the unintegrated gluon density in a proton which was obtained [38] from the CCFM 
equation. Of course, we expect that effects coming from the non-zero gluon virtualities for 
the associated bbH and tiH production are not very well prononced even at LHC energies. 
However, our study is important since it is planned to include the calculated off-shell matrix 
element g*g* QQH to the Monte-Carlo generator Cascade [39]. We will compare the 
results obtained in the fc^-factorization approach with the leading and next-to-leading order 
predictions of the coUinear factorization of QCD. 

The outline of our paper is following. In Section 2 we recall shortly the basic formulas of 
the /cy-factorization approach with a brief review of calculation steps. We will concentrate 
mainly on the g*g* —>■ QQH subprocess. The evaluation of qq — > QQH contribution is a 
rather straightforward and therefore will not discussed here (for the reader's convenience, 
we only collect the relevant formulas in Appendix). In Section 3 we present the numerical 
results of our calculations and a discussion. Section 4 contains our conclusions. 

2 Theoretical framework 

2.1 Kinematics 

We start from the kinematics (see Fig. 1). Let p^^"^ and p*^^^ be the four-momenta of the 
incoming protons and p be the four-momentum of the produced Higgs boson. The initial 
off-shell gluons have the four-momenta ki and k2 and the final quark and antiquark have the 
four- momenta pi and p2 and masses mg, respectively. In the proton-proton center-of-mass 
frame we can write 

p« = v^/2(l,0,0,l), = v^/2 (1,0, 0,-1), (1) 

where ^/s is the total energy of the process under consideration and we neglect the masses 
of the incoming protons. The initial gluon four-momenta in the high energy limit can be 
written as 

= xip(^) + kiT, k2 = X2P^^^ + k2T, (2) 
^The calculations [30, 32] were performed using the relevant on-mass shell matrix element. 
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where kiT and fc2T are their transverse four-momenta. It is important that k^^ = —kfj. ^ 
and k|j, = — ^2T 0- From the conservation laws we can easily obtain the following 
conditions: 

kiT + ^2T — Pit + P2T + Pt, 
X\\/s — miTC^^ + m2Te^^ + rriTe^ , (3) 

where y and niT are the rapidity and the transverse mass of the produced Higgs boson, piT 
and P2T are the transverse four- momenta of the final quark and antiquark, yi, y2, myr and 
m2T are their center-of-mass rapidities and transverse masses, i.e. m|j. = rriq -\- p|j.. 

2.2 Off-shell amplitude of the g*g* QQH subprocess 

There are eight Feynman diagrams (see Fig. 2) which describe the partonic subprocess 
9*9* ~^ QQH at aal order. Let ei and €2 be the initial off-shell gluon polarization vectors 
and a and h the relevant eight-fold color indices. Then the relevant matrix element can be 
presented as follows: 



KA 2-t V Pi-h + mi ki-p2 + m2 , . 

-Ma = / i.(Pi)*°7% 77 "^"".. g »(P2), (6) 

"^1 - (pi - "^f - {-p2 - pY 

KA 2-1 \ A, V pi-k2 + mi ^„ „ —p2—p + mi . , 
M4 = 9 u{pi)t'Y^.—2 — 7 TTT^^ 7^^M 2 / ^Hu{p2), 7 

KA 2 -/ N Ej- Pi +P + "^2 A;i-p2 + m2 ,„ „ . . 

= / fiO>0 «°7% + «'7% "(P.) . (9) 

T^2 - (Pl + P) "^2 - («^2 - P2)^ 



Mr = 5' «(Pi) YC^'^'i^. k2, -k, - k2) .''^:,, r'''t'^x 

[ki + k2Y 

^, -p2-p + mi 

Hu{p2), 



f^i - {-P2 -pY 

KA 2-( \ TT Pi+p + m2 
ml — (pi + pY 

xYc^-p{k,, k2, -k, - k2) ,. ,, r''t'uiP2). 

[ki + k2Y 



(10) 



(11) 



In the above expressions C'^^f{k,p,q) and H are related to the standard QCD three-gluon 
coupling and the if-fermion vertexes: 

C^^^ik^p, q) = gf^^ip - kf + g^^iq - pf + gP^{k - g)^ (12) 



4 



e mg^ (14) 
sin mz ' 

where is the Weinberg mixing angle and mz is the Z-boson mass. The summation on 
the initial off-shell gluon polarizations is carried out using the BFKL prescription [21-25]: 

Y.nh)^'^{h)-^^. (15) 

This formula converges to the usual expression X^e'^e^ — —g'^'' after azimuthal angle aver- 
aging in the A;t — > limit. The evaluation of the traces in (4) — (11) was done using the 
algebraic manipulation system Form [40]. We would like to mention here that the usual 
method of squaring of (4) — (11) results in enormously long output. This technical problem 
was solved by applying the method of orthogonal amplitudes [41]. 

The gauge invariance of the matrix element is a subject of special attention in the fcr- 
factorization approach. Strictly speaking, the diagrams shown in Fig. 2 are insufficient and 
have to be accompanied with the graphs involving direct gluon exchange between the protons 
(these protons are not shown in Fig. 2). These graphs are necessary to maintain the gauge 
invariance. However, they violate the factorization since they cannot be represented as a 
convolution of the gluon-gluon fusion matrix clement with unintegrated gluon density. The 
solution pointed out in [23, 24] refers to the fact that, within the particular gauge (15), 
the contribution from these unfactorizable diagrams vanish, and one has to only take into 
account the graphs depicted in Fig. 2. We have successfully tested the gauge invariance of 
the matrix element (4) — (11) numerically. 

2.3 Cross section for the QQH production 

According to the /cT-factorization theorem, the QQH production cross section via two 
off-shell gluon fusion can be written as a convolution 

a{pp ^ QQH) = I —fg{xi,klT, iJ,^)dklT^x 

X / ^f,{x2MT:l^')dklr^da{g*g* ^ QQH), 

where d'{g*g* — > QQH) is the partonic cross section, fg{x,'k'^,fi'^) is the unintegrated gluon 
distribution in a proton and 0i and 02 are the azimuthal angles of the incoming gluons. The 
multiparticle phase space Ild^Pi/2Eid^'^\j2p^'^~J2p'^^^) is parametrized in terms of transverse 
momenta, rapidities and azimuthal angles: 

^ = ^dp!^dy,^. (17) 
2Ei 2 ^ 2n ^ ^ 

Using the expressions (16) and (17) we obtain the master formula: 



xfg(xi,'k^T,IJ, )fg{x2,'k2T,l^ )dkiTdk2Tdp^TP2Tdydyidy2 



27T 271 27r 27r ' 
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where \Ai{g*g* QQH)\'^ is the off-mass shell matrix element squared and averaged over 
initial gluon polarizations and colors, ipi ^'^^ "^2 are the azimuthal angles of the final state 
quark and antiquark, respectively. We would like to point out again that \M.{g*g* — > 
QQH)\'^ strongly depends on the nonzero transverse momenta k^^^ and \<i^rp. If we average 
the expression (18) over 0i and 02 and take the limit k^^ and k^j. 0, then we recover 
the expression for the QQH production cross section in the collinear aal approximation. 

The multidimensional integration in (18) has been performed by means of the Monte 
Carlo technique, using the routine Vegas [42]. The full C-I--I- code is available from the 
authors upon requescl. 

3 Numerical results 

We now are in a position to present our results. According to (18), in the numerical 
calculations below we have used the CCFM-evolved unintegrated gluon density in a proton, 
namely set AO [38] . This set is widely discussed in the literatur^ and has been implemented 
in the Monte-Carlo generator Cascade [39]. As it often done for QQH production [17-20], 
we choose the renormalization and factorization scales to be equal: = fip = ^{f^Q + 
In order to investigate the scale dependence of our results we will vary the scale 
parameter ^ between 1/2 and 2 about the default value ^ = 1. As it was proposed in [38], 
for ^ = 1/2 and ^ = 2 we use the AO— and A0+ sets of unintegrated gluon densities, 
respectively. For completeness, we set to mi, = 4.75 GeV, rrit = 172 GeV, rriz = 91.1876 GeV, 
sin^ 9w = 0.23122 and use the LO formula for the coupling constant as(/i^) with nj = 4 
active quark flavours at Aqcd = 200 MeV, such that as(M|) = 0.1232. 

We begin the discussion by presenting our numerical results for the associated bbH and 
ttH total cross sections as a function of Higgs boson mass for the LHC energy, ^/s = 14 TeV. 
We consider 100 < ttih < 200 GeV since the production of a Higgs boson in association with 
a pair of beauty or top quarks at the LHC will play an important role only for relatively 
light Higgs bosons. The solid histograms in Fig. 3 correspond to the results obtained in the 
/cT-factorization approach of QCD with the CCFM-evolved gluon density. The theoretical 
uncertainties of these predictions are presented by upper and lower dashed histograms. The 
dash-dotted histograms represent results which were obtained in the standard (collinear) 
approximation of QCD at LO0. The contributions from the quark-antiquark annihilation 
mechanism, qq — > QQH, are shown by the dotted histograms. One can see that using of 
the fcy-factorization approach leads to some enhancement of the predicted bbH cross section 
at low rriH region (namely rriH < 150 GeV) in respect to the collinear LO QCD results. In 
the case of tiH production, the calculated cross sections in both approaches are very close 
to each other. It is because the large-x region, namely x ~ 0.1, is only covered here and 
therefore there is practically no effects connected with the small-x physics. It is important 
that we find our leading-order predictions fully consistent with the corresponding LO results 
presented in [17-20]. The small visible differense can easily come from the different quark 

^ lipatov@ theory, sinp . msu.ru 

^See, for example, review [29] for more information. 

^Numerically, we have used the standard GRV (LO) parametrizations [43] of collinear parton densities. 
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and gluon densitie^. 

In contrast with the NLO QCD results, the /Ej-factorization approach not reduces the 
strong scale dependence of corresponding LO QCD predictions, which has been pointed out 
in [17-20]. We conservatively estimate this theoretical uncertainty to be at most of order 
40 — 50% (see Fig. 3). Such scale dependence is significant, of course, and this fact indicates 
the necessarity of inclusion of the high-order corrections to the fc^-factorization formalism. 
So far the fcr-factorization is based on the leading-order BFKL or CCFM evolution equations. 
On the other hand, the kernel of BFKL equation has been calculated already at NLO [45], so 
that in the small-x regime the fc-r-factorization can be formulated at NLO accuracy also [46]. 
At moment, this problem is not solved and much more further efforts should be concentrated 
in this field. We only mention here that the leading-order fc^-factorization naturally includes 
the high-energy part of the NLO coUinear corrections. 

Our predictions for the transverse momentum and rapidity distributions of the Higgs 
boson as well as associated beauty or top quark Q are shown in Figs. 4 — 7. The distributions 
on the azimuthal angle distance between the H and Q as well as the Q and Q are shown in 
Figs. 8 and 9, respectively. These calculations were performed using uih = 120 GeV. The 
comparison of the /c-r-factorization approach to the collinear one shows the some broadening 
of the transverse momentum distributions due to extra transverse momentum of the colliding 
off-shell gluons. Also the fcy-factorization result shows a more homogeneous spread of the 
azimuthal angle A(j)^~'^ distance. At the same time, the cross sections calculated as a 
function of rapidities and as well as the azimuthal angle distributions da / dAcj)'^^^ 
show a similar behaviour, except for the overall normalization. 

To elaborate the difference between the ^T-factorization approach and the collinear ap- 
proximation of QCD, we investigate more exclusive observables, like the cross section differ- 
ential in the total transverse momentum of the QQH system, p^'^^ . In the usual collinear 
factorization of QCD the effect of intrinsic transverse momenta of the initial gluons can 
not be described until higher order corrections are taken into account. In the NLO QCD a 
non-zero pj.'^^ is generated by the emission of an additional gluon, while at LO it is always 
balanced to zero. In the /cr-factorization formalism, taking into account the non-vanishing 
initial gluon transverse momentum k^^ leads to the violation of back-to-back kinematics 
even at leading order. This effect is clearly illustrated in Fig. 10, where we plot the bbH 
and tiH cross sections as a function of pj!^^ . Note that only the off-shell gluon-gluon mech- 
anism, g*g* QQH, has been taken into account here. The relevant contribution from 
the quark-antiquark annihilation, qq — > QQH, is expected to be almost negligible for bbH 
production and probably can be sizeble for tiH on^. Keeping in mind that the NLO for this 
observable is the first non-trivial order, it would be useful to compare the NLO QCD and 
fc-r-factorization predictions in order to investigate the exact effect of high-order corrections 
in collinear factorization. 

In addition, we evaluate the fully exclusive cross section for bbH production by requiring 
that the transverse momentum of one or both final state beauty quarks be large than some 
p™* value. This corresponds to an experiment measuring the Higgs decay products along 
with one or two high pt beauty quark jets that are clearly separated from the beam. In 

^The LO parton densities from CTEQ5L set [44] have been used in [17-20]. 

^We do not consider here the problem of proper transverse momentum generation of initial state quarks. 
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Fig. 1 1 wc illustrate the dependence of these exclusive cross sections on the p™* parameter. 
Reducting the p™* from 50 GeV to zero approximately increases the relevant cross sections 
by a factors of about 10 and 100, respectively. In the coUinear factorization of QCD, if both 
beauty quarks are required to be produced with pt > 20 GeV, the NLO corrections reduce 
the LO predictions, and these corrections are positive if bcaiity quarks produced at small 
Pt [17]. However, one can see that predictions of the /cy-factorization approach overestimate 
the coUinear LO results in a wide p^* range. 

Finally, we would like to mention that our /cr-factorization calculations can be straight- 
forwardly generahzed to the case of scalar Higgs bosons of the MSSM by replacing the SM 
beauty and top quark Yukawa couplings with the corresponding MSSM ones. It is because 
the off-mass shell matrix element calculated above (see Section 2.2) is proportional to the 
beauty and top quark Yukawa couplings. In the MSSM, these couplings to the scalar Higgs 
bosons, gQQH, are given by a simple rescaling of SM couplings gggn [47], i.e. 

sin a ^ cos a 

9bbh0 — n 9bbh, QuhP — — — n 9tth 

cosp smp /..„s 

COS a ^ sina ^ ' 

QbbHO — Z dbbh-i gttH° — — ^ 9tth, 

COS p sm p 

where hP and H° are the lighter and heavier neutral scalars of MSSM, and a is the angle 
which diagonalizes the neutral scalar Higgs mass matrix. However, at the NLO level this 
rescaling is spoiled by one-loop diagrams in which the Higgs boson couples to a closed quark 
loop. We do not consider supersymmetric-QCD corrections in this paper. 



4 Conclusions 

We have studied the associated production of Higgs boson and beauty or top quark pair 
in hadronic collisions at the LHC conditions in the /cr-factorization approach of QCD. Our 
consideration is based on the amplitude of off-shell gluon-gluon fusion subprocess g*g* — > 
QQH . The corresponding off-shell matrix elements have been calculated for the first time. 
Sizcblc contributions from the qq QQH mechanism have been taken into account in the 
LO approximation of coUinear QCD. 

We have investigated the total and differential cross sections of hhH and ttH production. 
In the numerical calculations we have used the unintegrated gluon distributions obtained 
from the CCFM evolution equation. The comparisons with the leading and next-to-lcading 
order QCD predictions have been made. We demonstrate that the /c^-factorization ap- 
proach not reduces the strong scale dependence of coUinear LO QCD predictions, pointed 
out in [17-20]. This fact indicates the importance of the high-order correction within the 
fc^-factorization approach. These corrections should be developed and taken into account in 
the future applications. Finally, we show how our results can be generalized to the scalar 
Higgs sector of the MSSM. Our calculations is also important for Higgs boson searches where 
one or two high-pr beauty quarks are tagged in final state. 
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6 Appendix A 

Here we present the compact analytic expressions for the qq QQH subprocess. Let us 
define the four-momenta of the incoming and outgoing quark as fci, A;2, Pi and p2, respectively. 
The outgoing quarks have mass mq, i.e. p\—p\ — rrtq. In the formulas below we will neglect 
the masses of the incoming quarks. 

The contribution of the qq QQH subrocess to the total QQH cross section can be easily 
calculated using the master formula (18). One should only replace the unintegrated gluon 
densities /g( ) by the quark ones, perform the summation over initial quark flavours 

and take the coUinear hmit. The squared leading-order matrix elements \M.{qq — > QQH)\'^ 
summed over final polarization states and averaged over initial ones can be written as follows: 



\M{qq - QQH)\- = (^)^a^ 

72 sm 20w \^zf 



Fii F22 F12 + F21 



2 



{A.l) 



where 

Fii = 128 {piP2)iP2h){p2k2) - 64 {piP2){p2ki)ihk2) - 64 {piP2)iP2k2)ikik2)- 

64 {piP2){kik2)ml - 64 {piki){p2ki){p2k2) + 64 {piki){p2ki){kik2) + 
64 {pik{){p2k2f + 64 {piki){kik2)ml - 64 {pik2){p2ki){p2k2) + 
64 {pik2){p2k\f -\- 64 {pik2){p2k2){kik2) + 64 (pi/c2)(^i^2)"t.q- 

128 ip2ki)ip2k2)ml + 64 {kik2)m% + 64 {kik2fml, {A.2) 
F22 = 128 {piP2){piki){pik2) - 64 {piP2){piki){kik2) - 64 {piP2){pik2){kik2)- 
64 {piP2){kik2)'m^Q - 64 (piki) {pik2) {p2ki) - 64 {piki){pik2){p2k2)- 
128 {piki){pik2)ml + 64 {piki){p2ki){kik2) + 64 {piki f{p2k2) + 
64 {Pik2){p2k2){kik2) + 64 {pik2f{p2ki) + 64 {p2ki){kik2)ml+ 

64 {P2k2){kik2)ml + 64 {kik2)m% + 64 {kik2fml, {A2) 
F12 = F21 = -64 {piP2){piki){p2k2) - 32 {piP2)ipiki){kik2) - 64 {piP2)ipik2){p2ki)- 
32 {piP2){Pik2){kik2) - 32 {piP2){p2ki){kik2) - 32 {piP2){p2k2){kik2)- 
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64 {piP2){kik2)mQ + 64 (piP2)^(^i^2) - 32 ipiki)ipik2){p2ki)- 
32 {piki){pik2){p2k2) - 32 {piki){p2ki){p2k2) + 64 {piki){p2ki){kik2)+ 
64 {piki){p2k2)mQ + 32 {piki){p2k2f + 32 (piA;i)(/ciA;2)mQ+ 

32 {pikif{p2k2) - 32 {pik2){p2ki){p2k2) + 64 (piA;2)(p2A;i)mJ+ 
32 (piA:2)(p2fci)' + 64 (piA;2)(p2A:2)(A;iA;2) + 32 (piA;2)(A;iA;2)mJ+ 
32 {pik2f{p2ki) + 32 {P2ki){kik2)ml + 32 ip2k2){kik2)ml + 64 {k^k2fml, {AA) 
Ti^{p2ki) + {p2k2)-ikik2), T2^ipiki) + {pik2)-{kik2). {A.5) 
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Figure 1: Kinematics of the g*g* — > QQH process. 
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Figure 2: Feynman diagrams which describe the partonic subprocess g* -\-g* QQH at the 
leading order in as and a. 




14 



X 100 



10'^ I — ' — ' — ' — ' — ' — ' — ' ' — ' 
100 150 200 

m, (GeV) 

10^ E ' ' ' ' ' ' ' ' ' 3 




I , , , , , , , , , 1 

100 150 200 

m, (GeV) 

Figure 3: The total cross section of associated bbH (upper plot) and ttH (lower plot) pro- 
duction as a function of Higgs mass at i/s = 14 TeV. The solid histograms correspond to 
the results obtained in the fcr-factorization approach of QCD with the CCFM-evolved unin- 
tegrated gluon density (set AO). The upper and lower dashed histograms represent the scale 
variations of /cT-factorization predictions, as it was described in the text. The dash-dotted 
histograms represent results which were obtained in the standard (coUinear) approximation 
of QCD at LO. The contributions from the quark-antiquark annihilation mechanism (multi- 
plied by a factor of 100 in the case of bbH production) are shown by the dotted histograms. 
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Figure 4: The transverse momentum distributions da/dp^ of associated hhH (upper plot) 
and tiH (lower plot) production calculated at mn — 120 GeV and = 14 TeV. Notation 
of the histograms is the same as in Fig. 3. 
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Figure 5: The transverse momentum distributions da/dpf of associated hbH (upper plot) 
and tiH (lower plot) production calculated at mn — 120 GeV and = 14 TeV. Notation 
of the histograms is the same as in Fig. 3. 
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Figure 6: The rapidity distributions da/dy^ of associated bbH (upper plot) and ttH (lower 
plot) production calculated at uih — 120 GeV and ^/s — 14 TeV. Notation of the histograms 
is the same as in Fig. 3. 
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Figure 7: The rapidity distributions da j dy^ of associated \)bH (upper plot) and ttH (lower 
plot) production calculated at mu — 120 GeV and \fs — 14 TeV. Notation of the histograms 
is the same as in Fig. 3. 



19 





10^ 



_l I I I I I I I I I I L_ 







1 



A(^""' (rad) 



Figure 8: The azimuthal angle distributions da / dA(f)^~'^ of associated bbH (upper plot) and 
tiH (lower plot) production calculated at ttih — 120 GeV and ^/s — 14 TeV. Notation of 
the histograms is the same as in Fig. 3. 
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Figure 9: The azimuthal angle distributions da / d^(j)'^~^ of associated hbH (upper plot) and 
ttH (lower plot) production calculated at mu — 120 GeV and = 14 TeV. Notation of 
the histograms is the same as in Fig. 3. 
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Figure 10: The transverse momentum distributions da/dpf of associated hhH (upper 
plot) and tiH (lower plot) production calculated at mn = 120 GeV and = 14 TeV. The 
off-shell gluon-gluon fusion mechanism is only taken into account. 
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Figure 11: The cross sections for bbH production with one (upper histograms) or two (lower 
histograms) high-p^ beauty quarks as a function of the minimal 6-quark transverse momen- 
tum calculated at = 120 GcV and ^/s = 14 TcV. The solid and dashed histograms 
correspond to the results obtained in the /cy-factorization approach and in the coUinear LO 
approximation, respectively. 
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